Temperature Coefficient of Resistivity in Amorplious Semiconductors 
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By invoking the microscopic response method in conjunction with a reasonable set of approxima- 
tions, we obtain new expUcit expressions for the electrical conductivity and temperature coefficient of 
resistivity (TCR) in amorphous semiconductors, especially a-Si:H and a-Ge:H. The predicted TCR for 
n-doped a-Si:H and a-Ge:H is in agreement with experiments. The conductivity from the transitions 
from a localized state to an extended state (LE) is comparable to that from the transitions between 
two localized states (LL). This resolves a long-standing anomaly, a "kink" in the experimental logjQ c 
vs. T^^ curve. 

PACS numbers: 71.23.An, 71.38.Fp, 71.38.Ht. 



The temperature coefficient of resistivity (TCR) in dis- 
ordered systems is an important physical observable, 
is difficult to compute, and is of technological inter- 
est for microbolometer materials for thermal imaging 
applications||ll0]- Boltzmann or master equations are 
often used to calculate transport coefficients in crys- 
talline semiconductors and semi-metals. The low car- 
rier concentration in these materials results in a low ki- 
netic energy of carriers. Thus the Landau-Peierls condi- 
tion: the basic criterion for a kinetic approach [fi/r < 
max{Ef, fcgT}] cannot be satisfiedliSJ, where t is the 
time interval between two collisions (carrier with dis- 
order and /or phonon), and Ef is the Fermi energy of a 
material. Then, neither the elastic scattering by disor- 
der nor the inelastic scattering by a phonon has a well- 
defined transition probability per unit time. The situa- 
tion for amorphous semiconductors (AS) is even more 
difficult. Due to the strong electron-phonon (e-ph) in- 
teraction for localized states, any transition involving 
localized state(s) requires a reorganization of the vibra- 
tional configuration|j-|6|] . Energy conservation between 
the initial and final electronic states (s!] for these intrin- 
sic multi-phonon transitions is violated more seriously 
than for the single-phonon processes. Thus the kinetic 
method is unjustified for AS0]. 

In the kinetic approach^], it is often supposed that (i) 
electrical conduction is fulfilled by the transition from 
a localized state to another localized state (LL) and the 
transition from an extended state to another extended 
state (EE)[li 2|; (ii) the transition from a localized state 
to an extended state (LE) and the transition from an 
extended state to a localized state (EL) do not directly 
contribute to conduction; (iii) LE and EL transitions 
only maintain the non-equilibrium stationary distribu- 
tion of carriers between localized states and extended 
states during a conductiori_process. Although phonon- 
assisted delocalization|@, |^ and photon-excited tran- 
sient currentjioll have been considered intuitively, rig- 
orous expressions for LE and EL transition contributions 
to the conductivity are not yet available. 

Because the interaction between the external electro- 



magnetic field and an AS can be expressed with ad- 
ditional terms in the Hamiltonian, the transport coef- 
ficients can be expressed with transition arnplitudes in 
the Microscopic Response Method (MRM)jl,[l]. Thus 
the long time limit required in the kinetic approach^] 
is avoided. In addition, the MRM categorizes transport 
processes with diagrams computed to any given order 
of residual interactions t^. We have seen that even to 
zero order in the residual interactions, LE and EL tran- 
sitions contribute to conductivity Indeed, if one cal- 
culates the electrical conductivity of an AS from the full 
density matrix rather than its diagonal elements (mas- 
ter equation), one sees that LE and EL transitions con- 
tribute directly to electrical conduction. Since the MRM 
is equivalent to the density matrix method of Kubo lflU] , 
the two methods reach the same conclusion. 

Disorder scattering in EE transitions driven by field 
has been treated in the coherent potential approxima- 
tion. The conductivity from the EE transitions de- 
pends weakly on temp er ature fl2ll , and is the same order 
of magnitude as that from LL transitions above room 
temperaturelHH]. In this Letter we apply MRMjl,@] 
to derive the contributions to conductivity from LL, LE 
and EL transitions solely drive by an external field. Two 
examples, the conductivity and TCR of a-Si:H and a- 
Ge:H are described. 

An accurate conductivity calculation requires (i) the 
eigenvalues and eigenvectors of single-electron states 
and (ii) the eigenfrequencies and eigenvectors for each 
normal mode||6|]. To express the conductivity in terms 
of accessible material parameters, we approximate the 
vibrations of an AS by a continuous medium. Al- 
though translational invariance is destroyed in AS, the 
standing wave modes are still well-defined. Because 
most amorphous materials are isotropiclJ, [3] and only 
acoustic modes are important for the e-ph interaction in 
semiconductors lfisll , one can use a;^ = cA: for the vibra- 
tional spectrum, where a;^ is the angular frequency for 
the mode characterized by wave vector k. c is the av- 
erage speed of sound defined by 3 /c^ = + l/c^, 
where ct and c; are the speeds of transverse and Ion- 
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gitudinal sound waves, determined by the bulk and 
shear modulus. The Debye cutoff wave vector ku = 
(67r^M„)^/'' is determined by the number density «„ of 
atoms in an ASfl5]. Because the displacements of atoms 
satisfy a wave equation, the transformation matrix A be- 
tween the vibrational displacements ur at point R and 
the normal coordinates ©k characterized by wave vec- 
tor k is: 



A 



Rk 



(271) 
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(1) 



where V is the volume of a sample. For a-Si|0], c = 
6.21 X lO^m/s and fco = 1-44 x W^^m-'^. 

For simplicity, we take localized states to be spheri- 
cally symmetric llT3l . The difference between localized 
states is expressed by a single-parameter localization 
length ^^111]. We will use letter A with or without a nat- 
ural number subscript to denote a localized state. For a 
localized state A, denote Ra as the position vector of the 
center, the normalized wave function is 



-1/2^-3/2^^ 



(2) 



where r is the coordinate of the electron and local- 
ization length flj]. Following Mott, is determined 
by the eigenvalue E of localized state (^JaO: = 
{cZe^ / 4:TieQ£) {Ec — E)^^, where Z is the effective nu- 
clear charge of an atom core, £ is the static dielectric con- 
stant. Ec is the mobility edge of the conduction band, 
and c is a dimensionless constant. We will focus on 
n-doped material: transport in the conduction band, 
p-doped material may be treated analogously. For a- 
= 11.68, Z = 4 and the values of Ec are rather 
dispersedBSai: 0.2-2eV: we will takeO E, = 0.5eV. 
The most localized states are dangling bonds, they have 
the shortest possible = 2.35/2A (one half of a bond 
length). E = for a dangling bond, then c = 0.121. For 
many in the range of conduction band tail, the 

density of localized states (DOS) satisfies 



N(E) = («,„,/U)e-(^-^)/", 



(3) 



where U is the Urbach energy, nigc is the total num- 
ber of localized states per unit volume. For a-Si, U k 
50meV |l^[l9|| . the number density «/oc of localized con- 
duction states is 

^loc = 5/ (10.86A)''. The exponen- 
tial N(E) implies that most localized states in a-Si have 
a localization length in the range 6-12A. Denote n as the 
carrier concentration, the Fermi energy Ef of a weakly 
doped AS is: 



Ef = Ec + mn{n/2nioc), 



(4) 



When n < 2nigc, all occupied states are localized at 
T= OK. Ansatz (|3) only characterizes the band tail states. 



To describe dangling bonds, one can (i) introduce a rea- 
sonable DOS, e.g. a rectangle or a Gaussian; (ii) cor- 
respondingly modify Ef and energy zero-point for ex- 
tended states; (iii) add the contribution from the dan- 
gling bonds to Eqs. (|6l7b in the summation(s) over local- 
ized states. In this paper we ignore the small contribu- 
tion of dangling bonds to the conductivity. 

In an AS, an extended state is a packet of Bloch waves 
of its reference crystaljlillli, and is labeled by the wave 
vector of its principal Bloch wave, or more roughly by 
the momentum p of a plane wave 111, [l2h . For an AS, for 
which the reference crystal does not exist, a wave packet 
constructed from plane waves is still a reasonable ap- 
proximation for an extended state. We will use letter B 
with or without a natural number subscript to denote an 
extended state. Excepting EE transitions driven by an 
external field, we may approximate an extended state 
Xbi (r) by a plane wave with certain momentum p, and 
its eigenenergy is that of the plane wave: 



V 



^i/2gipT/ri 



Eb, = film, 



(5) 



where the zero-point of energy for extended states is at 
the mobility edge Ec. The attraction between an electron 
and an atom core may be approximated by a screened 
Coulomb potentialUl. For a-Si:H||l], we approximate 
its Thomas-Fermi wave vector by the value llSll for c-Si 
(?rjr = 1.7A~i. 

With the foregoing approximations, the velocity ma- 
trix elements in the expressions of conductivity can be 
computed [23]. One can also obtain the static displace- 
ments of the atoms in a localized state induced by the e- 
ph interaction and the reorganization energyfJ] for tran- 
sitions involving localized state(s)j3l, which are essen- 
tial input for the conductivity. 
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FIG. 1: (color online) logjQ cr vs. 1/T in two n-doped a- 
Si:H samples at a; = with carrier concentration 10^* and 
10^'cm^-'. The unit of a is ohm^^cm^^. 

We first calculate the conductivity from the LE tran- 
sitions (line 2b of Table 4 in Ref.||6|]). When ksT > fiU 
(the first peak of vibrational spectrum, hcJ = 232K for 
a-Si ll23l ), the two time integrals IbiA±, cari be approxi- 
mated by an asymptotic expansion |23|] , and the vibra- 
tional degrees of freedom are integrated out. For LE 
transitions, we first sum over final electronic states XIb 
for a fixed localized state A. It is convenient to use a 
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spherical coordinate system with R/^ as the origin and 
the wave vector direction k/|k| of the incident electro- 
magnetic wave as polar axis (z axis) i23ll . Because AS 
are isotropic, the angular part of the p integral can be 



carried out. One can show that a- 
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and 



•'yy 



a. Because (i) the center R^^ of (pA must be a 
neighbor of the observation point of current density, and 
(ii) the factors in the conductivity do not depend on R^^, 

^ ^f/_^^dE]V(E), where f = cZe'^ / {AneoeU) 



is the most probable localization length. For a-Si, 
11. 75 A, is quite close to the experimental value (2^ 
lOA. The conductivity from LE transitions is j23|] : 
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where 0^ = {MaLoJfiY'^O^^ and fiiOBA = Eb - Ea- 
^BA = \ YLa ^'^ctiOa)^ is the reorganization energy for 
transition (p^ — > One can show that Xba de- 

creases with £,A- 

Next we consider the conductivity from EL transitions 
driven by a field (line 6a of Table 5 in Ref.@]). Because 
the field-matter coupling is Hermitian, the conductiv- 
ity for EL transition driven by a field may be obtained 
from Eq.® by exchanging (p^ and Xb and noticing that 
Aab = '^BA ■ For the LE transition driven by transfer in- 
tegral and the EL transition driven by e-ph interaction, 
one does not have this symmetryjj,|6|]. 

To obtain the conductivity from LL transitions (line 
2a of Table 4 in Ref.@]), the velocity matrix elements 
iwj^'^^ and ^jj^^^ are computed with approximation ^ 
in a spherical coordinate system with as origin and 
R = Ra^ — Ry\ as polar axisj^]/ they exponentially de- 
cay with R = |R|. YIaAi can be carried out byiH first 
considering a fixed cp^ and scanning (p^^ at all possible R 
with different . The conductivity from LL transition 
driven by field becomes ll23ll 
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where Rc is the radius of physical infinitesimal volume 
elements fil]. ?^AiA = 2 En ^f^aC^a^ — ^a)^ is the reor- 
ganization energy for transition (pA — >• <pAi ■ To simplify 
notation, we used ^2 instead of ^^j, used ^1 instead of 

^A- 
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FIG. 2: (color online) TCR vs. T in two n-doped a-Si:H 



samples at a; = Hz with carrier concentration 10 
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10^'cm ^ , experimental values are taken from Refs. 1281.12' 



,-3 



and 



Eqs. l|6l7b express conductivity and its associated tem- 
perature dependence upon several material parameters: 
n, fiig^, U, Ec, e, qjF arid c. Rea as a function of tem- 
perature T is plotted in FiglT] for two n-doped a-Si:H 
samples||2. 127I1 with carrier concentration n = 10^^ and 
lO^^cm"^ [the unit of cr is ohm^^cm^^]. The LE conduc- 
tivity is the same order of magnitude as that from LL 
transitions, while the contribution from EL transitions is 
10^^ — 10^'' of that from LL transitions. Thus the cal- 
culated conductivity is the same order magnitude as the 
observed ones, cf. Fig. 8 of Ref. ll29ll . From Eqs. (|6l7b , one 
can easily compute the temperature coefficient of resis- 
tivity (TCR) X = jO^^jj = The correspond- 
ing TCR vs. T^^ is plotted in Fig|2l If there are several 
processes contributing to conductivity in a material, ac- 
cording to MRM, the total conductivity a of the material 
is tT" = Y^j (Tj, where <7j is the conductivity from the j 
process (LL, LE and EL ect.)[^. The overall TCR x re 



th 



iates to the TCR for each process by >c = cr ^J^j CTj) 



where = -ar^daj/dT. At 300-350K, the observed 
TCR is in range -0.02 to -0.08 for a n-doped a-Si:H with 



10i^cm-3j2i,|2l|3l- The calculated TCR from LL 



transitions is smaller than the experimental data, this is 
solid evidence that the contribution from LE transitions 
is important. 

In a-Si:H, a long-standing puzzle is (i) there is a kink 
in the observed log^Qa{T) vs. 1/T curve; and (ii) for 
very different doping concentrations, the kink tempera- 
tures collapse into a narrow range lsol] . We resolve these 
problems. The crossing temperature T* of (Tii{T) and 
'^le(T) is the key to understanding the kink. When 
T < T*, LL transitions are the main conduction mech- 
anism; for T > T*, LE transitions dominate conductiv- 
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FIG. 3: (color online)logjQ u vs. logjg « for n-doped a-Si:H and 
a-Ge:H. The unit of a is ohm^^cm^^, the unit of n is cm^^. 



ity. If one forced a single Arrhenius fit to the overall 
conductivity, the formal activation energy would be dif- 
ferent below and above T* . Thus one has a kink in the 
log]^gCr(r) vs. 1/T curve. The non-exponential behav- 
ior as indicated in Eqs.l l6l7b also has some role for the 
observed kink. For LL and LE transitions, we compute 
linear fits for the calculated logjQ <j(r) vs. T^^, for which 
the norms of the residuals are 0.03 (LL) and 0.15 (LE) for 



lO^^cm-^; 0.05 (LL) and 0.20 (LE) for lO^^cm"^. This is 
consistent with the deviation from linear relation in the 
measured mobility vs. 1/T curve, cf. Fig. 7.10 of 

From Fig |l(a)| and |l(b)[ we can see that T* de- 
creases with n: T* =294K for lO^^cm^^, r* =286K 
for lO^^cm^'^. This is consistent with the trend found 
in experiments: T* =400K for [PH3]/[SiH4]=10-6 and 
T* =333Kfor [PH3]/[SiH4]=10"2 cf pjg 3.1 of Ref.ii. 

Fig. |3(a)| and |3(b)| plot log^g a vs. log^g " T=273K 
and 300K. We see that at T=273K (300K), the contri- 
bution from LL transitions is larger (smaller) than that 
from LE transitions for carrier concentration from 10^^ 
to lO^^cm"^. For very different carrier concentrations 
(10^ times different), the kink temperatures fall near 273- 
300K, consistent with fact (ii). In Fig. |3(c) and |3(d)| we 



plot log^^Q (7 vs. log^o n for n-doped a-Ge:H at T=300K 
and 400K. T*s fall between 300-400K, are higher than 
those for a-Si:H. It agrees with the observations, com- 
pare Fig. 3.1 and 3.6 otfsd]. 

The compositional atomic orbital and / or their relative 
phases for the states close to Ep in the valence band (VB) 
are very different to those for states close to Ef in the 
conduction band (CB). Because U oc {[{Va - Vc)/Vc]^)av 
and Ec oc {{Vg — Vc)/Vc)av, where Va and Vc are the po- 
tential of an AS and its reference crystal[21], {)av denotes 
configurational and state average in AS|22]. Larger LI 



and Ec means stronger disorder, that implies smaller f 
or larger A, i.e. smaller u and larger k. 

Our approach is not restricted to one component sys- 
tems with weak electronic correlation. If one has a 
reasonable single-electron DOS in which the correla- 
tion between electrons is already taken into accoimt, it 
is not difficult to calculate the e-ph coupling in multi- 
component system, e.g. the e-ph interaction induced by 
the optical modes in VOi §3. The remaining procedure 
is exactly like here. 

In conclusion, the microscopic response method 
expresses transport coefficients with transition 
amplitudesj^ ^ rather than transition probability 
per unit time, which enables the method to be used 
with amorphous semiconductors for which the Landau- 
Peierls condition is violated||3|]. The conductivities from 
the three simplest transitions: LL, LE and EL transi- 
tions driven by field are expressed by several material 
parameters. The conductivity from LE transitions is as 
important as that from the LL transitions. The combi- 
nation is responsible for the kink in the experimental 
log^Q cr vs. T^^ curve. The LE transition is critical in 
determining the TCR. 

This paper provides new analytical form for a and x, 
suitable for amorphous semiconductors. A desirable ex- 
tension would be a full ab initio evaluation of all MRM 
diagrams using quantities from density functional the- 
ory. This complex task is underway. 
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